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ABSTRACT 

Aims. We present a forecast of dark energy constraints that could be obtained from a large sample of distances to Type 
la supernovae detected and measured from space. 

Methods. Wc simulate the supernova events as they would be observed by a EUCLID-like telescope with its two imagers, 
assuming those would be equipped with 4 visible and 3 near infrared swappable filters. We account for known systematic 
uncertainties affecting the cosmological constraints, including those arising through the training of the supernova model 
used to fit the supernovae light curves. 

Results. Using conservative assumptions and Planck priors, we find that a 18 month survey would yield constraints 
on the dark energy equation of state comparable to the cosmic shear approach in EUCLID: a variable two-parameter 
equation of state can be constrained to ~ 0.03 at z ~ 0.3. These constraints are derived from distances to about 13,000 
supernovae out to z = 1.5, observed in two cones of 10 and 50 deg 2 . These constraints do not require measuring a 
nearby supernova sample from the ground. 

Conclusions. Provided swappable filters can be accommodated on EUCLID, distances to supernovae can be measured 
from space and contribute to obtain the most precise constraints on dark energy properties. 

Key words, cosmology: cosmological parameters - cosmology:dark energy 



1. Introduction 

About ten years ago, distances to abo ut 50 Type la su- 
pernovae (SNe l a ) ena bled two teams (|Riess et al.l Il998t 
iPerlmutter et al.l I1999D to independently constrain the 
kinematics of the expansion of the universe and present the 
first evidence for acceleration at late times (typically later 
than z ~ 0.5). This acceleration was ascribed to a mys- 
terious component baptised dark energy. The mere cause 
of this late acceleration remains unknown, and the expan- 
sion history is not well enough constrained yet to uniquely 
characterise the phenomenon. There is no known candi- 
date to incarnate this dark energy, and working groups have 
been constituted to review possible observational strategies 
that could significantly increase our knowledge of it. Two 
work ing group reports (|Albrecht et al.ll2006t iPeacock et al.l 
2006) have identified a set of observational approaches to 
constrain dark energy, and outlined generic large and dim- 
cult projects, mainly (but non only) space-based, that could 
significantly improve our knowledge of dark energy. 

Four main dark energy probes were identified: the cos- 
mic shear correlations as a function of angle and redshift, 
the measurement of the acoustic peak in the galaxy correla- 
tion function, the measurement of distances to Type la su- 
pernovae, and cluster counts. All these measurements have 
to be carried out as a function of redshift in order to con- 
strain dark energy. Forecasts and merits of these methods 
were discussed in both reports which stress that crossing- 
methods is mandatory, since the anticipated measurements 
face partly unknown systematic uncertainties. Regarding 
the merits of the considered probes, a short summary of 
the findings is that cosmic shear correlations are the most 



promising and most demanding approach, baryon acoustic 
oscillations (BAO) is the simplest (in terms of analysis), 
and distances to supernovae the most mature. One strong 
incentive to use multiple probes, beyond redundancy, is 
that General Relativity predicts a specific relation between 
the expansion history and the growth of structures that can 
be uniquely tested by comparing expansion history from 
supernovae and BAO, to growth rate from cosmic shear, 
cluster counts and redshift distortions. (The measurement 
of redshift distortions can be extracted from the same data 
as BAO.) 

The Hubble diagram of Type la supernovae remains 
today the best dark energy probe, at least if one focuses 
on results rather than forecasts. Ambitious space-based su- 
pernova programs have been imagined and presented right 
after the discovery of accelerated expansiorQ that would 
observe 0(1000) supernovae from space. More recently, a 
small visible-imaging mission concept DUNE has been con- 
sidering the obser vation of about 10,00 supernovae up to 
z = 1 from space (|Refreeier et al.ll2006l ). This mission con- 
cept was later extended to the near inf rared (NIR) an d 
proposed to ESA under the same name (iRefregieri 12001 . 
New space missions are currently being developed in or- 
der to constrain dark energy, both in north America (the 
JDEM mission, now called WFIRST) and in Europe (the 
EUCLID mission, which incorporates the second DUNE 
concept). We propose here a space-based supernova sur- 
vey that could plausibly be implemented on the EUCLID 
project, provided a filter wheel is accommodated into the 
visible camera. For the first DUNE concept, suppressing the 



see e.g. http://snap.lbl.gov/ 
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filter wheel had been considered and had a negligible impact 
on the mission cost and complexity. The supernova survey 
forecast we present here relies on the exp erience gained in 
analysing the grou nd-based SNLS survey (|Guv et al.ll2010t 
iConlev et alll2010h . 

We will hrst describe the instrument suite we simulate. 
We follow in fj3] with the proposed supernova survey. The 
proposed methodology for the analysis and the associated 
Fisher matrix are the subjects of 21 Our results are pre- 
sented in fJSJ We discuss in ^5. 21 how the key issues related 
to distance biases might be addressed. We summarise in §5\ 



2. Instrument suite 

In its current design, EUCLID is equipped with a visible 
and a near infrared (NIR) imager (|Laureiisll2009t ). (It also 
embarks a slitless NIR spectrograph, and the SN survey 
we discuss here does not rely on it). We assume that both 
imagers observe simultaneously the same part of the sky 
through a dichroic, but if they were instead observing con- 
tiguous patches, it would marginally affect the efficiency of 
the proposed surveys. 



2.1. General features 

We assume the following figures and features : 

— The primary mirror diameter is 1.2m, with a central 
occultation of 0.5m. 

— Both the visible and NIR imagers cover 0.5 deg 2 . The 
survey efficiency driver is the field of view of the NIR 
imager. 

— The pixel sizes are 0.1" for the visible channel and 0.3" 
for the NIR channel. 

— We assume a quantum efficiency of 0.8 for the NIR de- 
tectors. The CCD sensors for the visible channel are 
assumed to have a quantum efficiency peaking at 0.93 
at 600 nm, and falling to 0.5 at 900 nm. The red sensi- 
tivity is probably conservative. 

— We assume that the NIR imager has a 20 electrons read- 
out noise and a 0.1 el./s/pix dark current. Both numbers 
are not particularly optimistic. 

— The visible imager has a 4 electrons read-out noise and 
a 0.002 el./s/pix dark current. Both have a negligible 
impact on broad-band photometry of faint sources. 

— At variance with the current EUCLID design, we as- 
sume that both imagers are equipped with filters, with 
a filter changing device such as a filter wheel. 

— The filters transmit at most 90% of the light, and the 
optics transmission (excluding filter and CCD) is about 
80% at 1.2 fim. 

— The image quality (IQ) is due to two components : one 
independent of wavelength equal to 0.17"(FWHM) and 
a diffraction contribution for a 1.2 m mirror. Within 
this model, the FWHM IQ increases from 0.2" at 0.5 
/xm to 0.35" at 1.5 /xm. The contribution of pixel spatial 
sampling will be considered later. 

— The light background is taken from iLeinert et al.l 
(1998). We assume that the survey takes place at 60 
degrees of ecliptic latitude, which corresponds roughly 
to a 20 % increase of background light with respect to 
the ecliptic pole. 
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Fig. 1. Overall transmission of the bands of both imaging 
systems as they are simulated for this study. The EUCLID 
baseline set consists of 3 NIR filters (as assumed here) and 
a single "r+i+z" visible filter. 



2.2. Filter set 

Cosmological constraints using supernovae derive from 
comparing measurements across redshifts. In order to make 
the comparison as robust as possible, we propose a filter set 
logarithmic in wavelength, with 7 bands: 4 bands in the vis- 
ible channel, c orresponding roughly to g,r,i and z from the 
SDSS system (jFukueita et al.l ll996). and 3 bands on the 
NIR channel roughly matching y, J and H. The NIR filters 
are already included in the EUCLID concept, but the visi- 
ble filters are not. The whole filter set covers the [450, 1660] 
nm interval, with the dichroic split at 950 nm. The overall 
transmission of the simulated bands is displayed in Fig. [TJ 
The sky brightness we assume in our simulations and the 
collecting power of the instrument (provided as magnitude 
zero-points) can be found in Table [TJ 



band 


\{A) 


sky (AB/arcsec 2 ) 


zp (c /s) 


10a- (600s, ps) 


g 


4965 


23.19 


24.27 


25.1 


r 


5958 


22.89 


24.32 


25.1 


i 


7162 


22.68 


24.18 


24.9 


z 


8574 


22.60 


23.72 


24.4 


y 


10476 


22.47 


24.32 


24.4 


J 


12549 


22.44 


24.37 


24.4 


H 


15148 


22.31 


24.41 


24.4 



Table 1. Characteristics of the bands considered in 
this work: central wavelength, sky brightness (in AB 
magnitudes/arcsec 2 ), zero-points (for AB magnitudes and 
fluxes in e~/s), and AB magnitude limits for point sources 
detected at S/N=10 in 600 s using PSF photometry, ac- 
counting for dark current, read-out noise and pixel sam- 
pling. Note that the sensitivities provided in lLaureiisI ([2009) 
differ substantially because they apply to objects measured 
using aperture photometry, with a higher sky background 
because the EUCLID wide survey cannot be confined to 
high ecliptic latitudes. 
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3. Characteristics of the supernova survey 

Supernovae are transient events and the primary distance 
indicator is the amplitude of the light curve(s), that we as- 
sume to be measured via broad-band imaging. In order to 
properly sample supernova light curves, one needs to re- 
peatedly image the same pointings, in as many bands as 
possible. We require that the rest frame wavelength region 
from 4000 to 6800 A is measured for all events, and sam- 
pled by (at least) three filters: this roughly corresponds 
to measuring the B,V and R (rest frame) standard bands. 
Current distance estimators only require two bands (they 
can of course accommodate more bands). Requiring three 
bands will either allow one to use more elaborate distance 
estimators, or to use the built-in redundancy to compare 
derived distances. Requiring that the same rest frame wave- 
lengths are measured at all redshifts make the comparison 
of supernovae brightnesses as independent as possible from 
a supernova model. Although we require three bands to 
measure a distance, the observing strategy proposed below 
provides us with more than three bands at all redshifts. The 
duration of the survey is arbitrarily chosen as 1.5 year, i.e. 
about one third of a five-year mission. This range of survey 
time provides O(10 4 ) events, i.e. an order of magnitude ove r 
current cosmological samples (e.g. lAmanullah et alJl20lj) . 
Alterations of the survey duration will be considered in 



3.1. Photometric precision requirements 

The depth of the observations is driven by the maximum 
target redshift, and considering three requirements: (1) the 
measurements should be accurate enough so that the dis- 
tance estimate has a measurement uncertainty that re- 
mains significantly smaller than the observed scatter of the 
Hubble diag ram (namely around 0.15 mag for modern sur- 
IGuv et alJ2007HKessler et alJ200fllcTuv et al.1 



veys, seee.g.iv .u . . .... _ . . .. ... — ....... 

12010: iConlev et all 12010ft : (2) the same supernova param- 
eter space should be observable at all redshifts in order 
to avoid the shortcomings of Malmquist bias corrections; 
(3) the amplitude of the three bands roughly matching 
rest frame B,V,R bands should be measured at an accu- 
racy better or comparable to the measured " c olour smear- 
ing", namely about 0.025 mag, see lGuv et all (|2007l [2010h . 
Colour smearing refers to the spread of colour-colour rela- 
tions of SNe la: it is the scatter of single-band light curve 
amplitudes one has to add to measurement uncertainties in 
order to properly describe the observed spread of colour- 
colour relations of SNe la. One can for example model 
the rest frame V-R colour from B-V. The scatter of this 
colour-colour relation is properly described by assuming 
that the B,V and R peak magnitudes scatter independently 
by <j c ~ 0.025 mag around a two-parameter model. One 
can est imate the colour smearing as a function of wave- 
length (jGuv et alJ[201Q[ ). and find that the B,V,R region 
is less scattered than bluer bands. The colour smearing 
contributes to the Hubble diagram scatter, but does not 
account for the entirety of the measured ~0.15 mag r.m.s. 

In practice, the requirement (2) about Malmquist bias 
indicates that events one magnitude fainter than the aver- 
age at the highest redshift should be easily detected. We 
find that requirement (3) is the most demanding, and the 



survey setup we present later fulfils these three require- 
ments. 



3.2. Redshifts 

Supernova redshifts are obviously needed to assemble a 
Hubble diagram. Obtaining SNe la spectra from the ground 
is technically feasible up to z ~ 1. At higher redshifts, the 
required exposure times increase very rapidly with redshift 
because the SNe la flux decreases rapidly in the near UV 
(bluer than ~ 3600 A), and the atmosphere glow increases 
rapidly towards the red. A large spectroscopic followup of 
supernovae at z > 1 is hence out of reach of current ground- 
based instruments, and even if adaptive optics, OH sup- 
pression and larger telescopes will certainly help very sig- 
nificantly, massive statistics will likely remain out of reach 
during at least the next decade. 

Since acquiring a spectrum for each supernova was a 
core goal of the proposed SNAP conceptQ its baseline 
design included a high-throughput low-resolution spectro- 
graph. In this concept, the supernova spectroscopy require- 
ments drove the mirror size and limited the statistics to 
0(2000) supernovae, because of the long exposures required 
for spectroscopy of distant supernovae. 

We believe that obtaining supernova spectra one at a 
time with a small field instrument is not a realistic goal 
for O(10 4 ) events, especially if aiming at redshifts beyond 
unity. Wide field space-based slitless spectroscopy may be 
considered, but the S/N ratio of slitless spectra is natu- 
rally poor, because every pixel integrates the sky back- 
ground spectrum in the whole bandwidth: spectra at a 
modest magnitude H ~ 24 seem out of reach of EUCLID 
slit-less spectrograph. Space-based spectroscopy with syn- 
thetic slits (as originally proposed for the SPACE project, 
later merged with DUNE into EUCLID) is significantly 
more sensitive, and for a supernova program on a joint 
imaging-spectroscopy-with-slits mission (such as the orig- 
inal EUCLID concept), one should obviously consider ob- 
taining supernovae spectra in parallel with repeated imag- 
ing. 

So, for O(10 4 ) supernovae, we should consider the case 
where supernova spectra cannot be acquired for all events. 
Since the Hubble diagram obviously requires redshifts, we 
consider the following alternatives: 

— Acquiring host galaxy spectra "after the fact" using 
wide field spectroscopy. The target density would be be- 
low 500 deg -2 and hence perfectly suited to the multi- 
fiber spectrography being considered today (WFMOS 
and BigBoss are typical examples). For the redshift 
range and area we will consider below, this would typ- 
ically require a few hundred hours of integration on a 
4-m class telescope, which seems acceptable, thanks to 
the multiplex factor. 

— Relying on photometric redshifts of host galaxies, using 
both visible and NIR bands. This would cause a loss 
of accuracy of cosmological parameters. However, one 
would still have to collect a sample of spectroscopic red- 
shifts of faint galaxies in order to train the photometric 
redshift algorithms. 

— Relying on photometric redshifts of supernova. These 
are now known to be more accurate than photometric 



The expression was introduced in Kcsslcr ct al. (200 
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redsh i fts of host galaxies (iPalanaue-Delabrouille et al.l 
120101 : iKessler et al.l l2010al) , thanks to the homogene- 
ity of the events. However, this approach weakens the 
identification step since the redshift will minimise the 
difference between measurements and expectations for 
a Type la supernova. It also introduces correlated un- 
certainties between distance and redshift which would 
require a careful study. 

We assume in what follows that supernovae have high 
quality redshifts, i.e. that spectroscopic redshifts of host 
galaxies are acquired. From here on, we will allow for a 
conservative 20 % loss in supernovae statistics due to this 
process. 

3.3. Why space? 

Two aspects favour a space based supernova survey: first, 
ground-based surveys suffer from variations of atmospheric 
transparency and image quality across observing epochs. 
This causes photometric uncertainty floors which can be 
overcome from space. Second, photometry of faint sources 
in the NIR from the ground is notoriously difficult and re- 
quires extreme exposure times on 8-meter class telescopes, 
because of the bright atmospheric glow. The "1 /an barrier" 
practically limits supernova ground-based surveys to z = 1 
or significantly below, if one imposes that supernovae at all 
redshifts are measured in the same rest frame wavelength 
range. Space offers NIR photometry with sensitivities simi- 
lar to the visible range, and hence opens supernova surveys 
to the z > 1 range. Accessing NIR bands also improves su- 
pernova surveys in two respects: it allows one to measure 
objects on a large wavelength lever arm, which is manda- 
tory to characterise their colour variations. It also allows 
one to measure distances to supernovae using redder bands 
than allowed from the ground (we will consider rest frame 
I band in H3.6I) . which limits the effects of colour varia- 
tions among the sample. Finally, accumulating deep NIR 
photometry of galaxies improves galaxy typing and photo- 
metric redshifts (over visible-only measurements); the lat- 
ter may be useful even for supernova cosmology if obtaining 
spectroscopic redshifts of the supernova hosts turns out to 
be not practical. 

3.4. Supernova surveys cadence and coverage 

The imaging survey should be run in "rolling search mode" 
where a given patch of the sky is observed repeatedly, in or- 
der to discover and measure variable objects. The required 
wavelength range to cover is bound to the redshift range 
aimed at: in order to derive distances as independent as pos- 
sible of any supernova model, the same rest frame spectral 
region should be used to derive distances at all redshifts. 

Monitoring a single cone at the depth required at the 
survey highest redshift delivers a supernova redshift dis- 
tribution where moderate and low redshifts are essentially 
missing. We hence propose a two-cone approach: a deep 
survey of 10 deg 2 up to z — 1.55, and a wide survey of 
50 deg 2 up to z = 1.05. We consider different area ratio 
in t j5.ll The deep survey typically requires integrations 4 
times longer than the wide. Note that both surveys are on 
purpose volume limited in order to avoid the shortcomings 
of selection bias corrections. 



We simulate a survey duration of 1.5 year of calendar 
time, with the wide and deep survey observations inter- 
leaved. The footprint of the deep survey is assumed to be 
outside the wide, so that its low redshift sample (where 
statistics is precious) adds up to that of the wide. We 
draw the simulated samples from a measured SNela vo- 
lum ic rate as a functio n of redshift (jRipochel (|2008h . see 
also lPerrett et al1l2010t ). which may be parametrised as 

R(z) = 1.53 10~ 4 [(1 + z)/1.5] 2 ' 14 ftfo Mpc- 3 yr" 1 

where years should be understood in the rest frame. Since 
these measurements stop around z = 1 and rates at larger 
redshifts are highly uncertain, we assume that the volumic 
rate remain s constant above z = 1 (to z = 1.5). The rates 
proposed in iMannucci et al.l (|2007f ) (accounting for events 
"lost to extinction") yield a statistic (to z = 1.5) ~ 25 % 
larger than what we simulate, with a similar redshift dis- 
tribution. Our simulation accounts for "side effects" by re- 
quiring that all epochs corresponding to rest frame phases 
from -15 to +30 days from maximum light are measured. 
As a consequence the event statistics increases if the sur- 
vey monitors smaller areas over a longer period, with a 
constant total observing time. The accuracies we discuss 
later ignore 20% of the events, in order to allow for losses 
in the measurement process (failures to obtain redshifts for 
example). Accounting for these losses, the deep and wide 
surveys deliver about 4000 and 9000 events respectively. 
The simulated statistics are provided as a function of red- 
shift in Table [2j Events at higher redshifts can be detected 
but they fail the quality cuts for deriving distances. 





^max 


wide 


deep 


0.15 


0.25 


120 


24 


0.25 


0.35 


263 


52 


0.35 


0.45 


454 


90 


0.45 


0.55 


685 


137 


0.55 


0.65 


945 


189 


0.65 


0.75 


1226 


245 


0.75 


0.85 


1520 


304 


0.85 


0.95 


1819 


363 


0.95 


1.05 


2090 


418 


1.05 


1.15 




435 


1.15 


1.25 




440 


1.25 


1.35 




442 


1.35 


1.45 




439 


1.45 


1.55 




434 


total 


9122 


4012 



Table 2. Expected number of SNe la events in 1.5 year for 
the wide and deep surveys respectively (accounting for a 
spectroscopic efficiency of 0.8). Events beyond a redshift of 
1.05 in the wide survey are ignored because the quality of 
their measurements is below our requirements. 



We settle for a cadence of 5 observer days, but the pho- 
tometric accuracies that matter to derive distances do not 
depend on this choice at first order, since those mainly 
depend on the overall integrated light. However, a signif- 
icantly coarser sampling might compromise photometric 
identification. The integration times per filter at each epoch 
are provided in Table |3] and were chosen in order to pro- 
vide a measurement of light curve amplitudes to a precision 
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of 2.5 % r.m.s on average at the highest redshifts of each 
survey. 



Filter 


wide 


deep 


g 


400 


1600 


r 


400 


1600 


i 


500 


1600 


z 


700 


2400 


y 


650 


1800 


J 


650 


2400 


H 


650 


3000 



Table 3. Integration times per visit (in seconds) proposed 
for the wide and deep supernova surveys with a 5 days 
cadence. The wide and deep surveys have to cover respec- 
tively 100 and 20 pointings. The CCD type we simulate has 
a modest quantum efficiency in the z band, which drives ex- 
posure times up in this band. 



In 5 days of calendar time, the total exposure times 
amount to 56 and 40 hours for the wide and deep surveys 
respectively, which corresponds to 80 % of the 120 hours 
available. 

Both surveys are conducted in 7 bands, covering g to 
H (450 < A < 1700 nm). One may argue that the bluest 
bands of the deep survey and the reddest bands of the wide 
survey are not strictly needed to measure supernovae dis- 
tances, but since visible and NIR observations are assumed 
to happen in parallel, dropping blue visible bands for the 
deep survey, or red NIR bands for the wide survey does not 
save any observing time. 

Note that we have not discussed the collection of a 
nearby sample (typically 0(1000) events at 0.03 < z < 0.1). 
As a baseline, and at variance with the SNAP project, we 
stick to a self-contained imaging survey, in order to realis- 
tically limit cross-calibration and detection bias issues. 

3.5. Identification of SNe la and contamination 

SNe la exhibit reproducible res t frame colours (< j(B — 
V) ~ 0.1, see e.g. Fig. 8 of lAstier et all I2006D . and 
even more reproducible colour relations. For example, the 
rest frame U-band amplitude can be predicted to bette r 
than ~ 0.04 from B- and V-bands (j Astier et al.l l2006h . 
SNe la not only occupy narrow subspaces of multicolour 
spaces, but also exhibit very reproducible light curve shapes 
that permit to discriminate against most of the core- 
collapse events (|Poznanski et al.l l2002t iJohnson fe Crottsl 
2006; iRodnev fc Tonrvl 12009b "studies of the photometric 
selection of SNe la have been conducted on the S NLS data 
and their prelimin ary conclusions are encouraging (|Ripochd 
2007; Bazin 2008). These studies rely on host galaxy photo- 
metric redshift and most of their identification failures are 
due to wrong assumed redshifts. The supernova survey we 
are considering here would be in a more favourable situa- 
tion than these studies: it measures 7 bands (whilst SNLS 
has at most 4), and we assume that host galaxy spectro- 
scopic redshifts will be available (whilst SNLS studies used 
host galaxy photometric redshifts). 

In order to estimate the contamination by core-collapse 
supernova in a sample selected in colour-colour subspaces, 
we would need a large enough sample of multi-band mea- 



surements of such events. These should soon be available, 
thanks at least to the SDSS supernova survey, the Lick 
Observatory Supernova Search, and the Palomar Transient 
Factory, but they are not available yet. Therefore, in or- 
der to bound the impact of core-collapse contamination on 
the SNe la distance-redshift relation, we resort to studying 
how clipping around the Hubble l i ne rej ects other super- 
nova types, following IConlev et al.l (|2010h . 

Core collapse supernovae are classified in Type lb and 
Ic, and Type II. Type lb and Ic events are often merged 
into a "Ibc" type (see e.g. iRichardson et al.ll2002t iLi et al.l 
120 101). Type II eve nts are about 3 times more frequent than 
Ibc (|Li et aDl2010l Fig. 9), but two thirds of those are Type 
Il-plateau (II-p) which are easily identified from their very 
flat light curves. Other Type II events represent about the 
same rate as Ibc, but their light curves rise in a few days, 
whereas SNe la rise in more than 15 days. So, Type II 
events might add a small contribution to Ibc interlopers, 
and we will now concentrate on evaluating the impact of a 
Ibc contamination in the Hubble diagram. 

We model the Ibc population absolute magnitude dis- 
tribution as a Gaussian offset by A;, c from the la pop- 
ulation with r.m.s cr£, c . We expect the rate of Ibc to be 
proportional to the s tar fo rmation rate, that we take from 
Hop kins fc Beacom (2006), and the amount of Ibc events 
follows from fb c (z — 0), the ratio of the Ibc rate to the la 
rate at z — 0. The adopted la rate was presented in ^3.41 

We simulate a mix of la and Ibc events, fit the Hubble 
diagram with a sixth degree polynomial, and iteratively clip 
events beyond 3 <j from the fit and refit, until no event 
is clipped. We simulate la events with a Gaussian scat- 
ter 0.15 mag around the Hubble line. For Ibc events, it 
is unclear how applying to them the brighter-slower and 
brighter-bluer corrections for SNe la will affect the absolute 
magnitude distribution. We might guess that some part of 
their brightness scatter is due to extinction in their host 
galaxy, and that brighter-bluer corrections would narrow 
their distance modulus distribution. We will consider below 
two estimates of the Ibc brightness scatte r, one as observed 
and o ne significantly lower. Note that IRichardson et all 
(2002) estimate intrinsic magnitude scatter corrected for 
host galaxy extinction which differ little from raw esti- 
mates. The contamination depends on the bright end of 
the Ibc luminosit y func tion which is not well constrained: 
IRichardson et al.l (|2002h propose a distrib ution extending 
well beyond the la average brightness, while lLi~eT al. (2010) 
brightest Ibc event is 0.5 mag fainter than the average la. 
Both results are however compatible at the ~10 % CL, 
given the modest statistics involvecQ- We propose three sce- 
narios for the Ibc population, each with two values for er;, c , 
detailed in Table HJ 

— (R02) the Ibc population in IRichardson et al.l (|2002h 
amounts to about 16 % of the SNe la, is fainter by 1.4 
mag than S Ne las and has a r.m.s. s catter of 1.4 mag. 

— (R02 bright) IRichardson et all (|2002h see a mild indica- 
tion of a bright component of Ibc, representing about 5 
% of the SNe l a. 

— (L10) ILi et al.l (|2010h have a much more complete sur- 
vey and find a volumetric rate of Ibc which is about 80 

4 Assumin g the disagreem ent is real, it might be due to the 
fact that the lLi et alJ (|20 lOT l sample comes from a search target- 
ing nearby galaxies and could be missing events preferentially 
occurring in dwarf galaxies. 
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% of the la rate. Their average Ibc is fainter by 2.4 mag 
than their average la with an r.m.s of 1.2 mag. as we 
just noted, the magnitude distribution of the measured 
sample does not contain Ibc events under the la peak. 

Contamination affects cosmology by biasing the average 
distance modulus. However, a redshift independent bias has 
no effect on cosmology, and we report in Table 0] the slope 
of the bias as a function of redshift dSfi/dz, which fairly 
describes most of the effect. 



Scenario 




fbc 


Obc 


d Sn/dz 


N bc /N Ia 


R02 


1.4 


0.16 


1.4 
1.0 


0.0014 
0.0034 


5.5% 
4.9% 


R02 bright 


-0.8 


0.05 


0.5 
0.33 


-0.0047 
-0.0047 


2.6% 
1.6% 


L10 


2.4 


0.8 


1.2 
1.0 


0.0086 
0.0058 


8.4% 
4.3% 



Table 4. Various hypotheses about the average magnitude 
shift of Ibc's with respect to SNe la, the r.m.s of their mag- 
nitude distribution, the amount of Ibc at z = in units of 
the la rate, the resulting slope bias of the Hubble diagram 
d Sfi/dz, and the average fraction of contaminants. The im- 
pact on cosmology is small under all these hypotheses. 



We find that under the three scenarios, the effect on 
cosmology is small. Even in our "L10" scenario, where the 
contamination is suspiciously large compared to the frac- 
tion of Ibc identified by high redshift SNe la surveys, the 
effect on cosmology remains below the systematic uncer- 
tainty <j{cm) — 0.01 (defined later in Eq. [SJ that we will 
consider as a baseline. The clipping process has a negligible 
impact on the SNe la statistics. A rough estimate of the 
size of the effect we find can be readily computed, for e.g. 
the first line of the table: the average offset of Ibc within 
the la ±3<T window is 0.047, the fraction of the Ibc popu- 
lation within the same window amounts to ^20%, the Ibc 
to la number ratio at z = is assumed to be 0.16, and it 
increases with redshift as ~ (1 + z). The evolution of the 
Ibc induced distance bias reads 0.047 * 0.2 * 0.16 = 0.0015. 

One might argue that a Gaussian distribution is inade- 
quate to describe the tails of the Ibc distribution, but more 
populated tails also exhibit shallower slopes, and both al- 
terations have a tendency to cancel each other. Note that 
our estimates ignore the potential rejection from light curve 
shapes and colours of Ibc events, which is certainly a con- 
servative assumption. 

Our r esults might look at odds with other attempts. 
iHom eicr (2005) follows a similar approach and finds sizable 
effects on cosmology, and we attribu te the differen c e to the 
absence of clipping. More recently, iKessler et al.l (|2010bj ) 
proposed a supernova photometric identification challenge. 
The provided simulated sample contains a large fraction of 
events for which the phase coverage and the signal to noise 
ratio are inadequate for a distance measurement, would 
they be genuine Type la. In our simulation, such events 
are excluded a priori and do not count as identification 
failures. We also benefit in our simulation from assuming 
spectroscopic redshifts as this reduces the Hubble diagram 
scatter. Regarding the iden tification usin g colour s, Type 
la events were generated in IKessler et al.l (|2010b[ ) with a 
colour smearing of 0.1 mag, which we now know to be 



about 4 times too large, and leads to overestimating mis- 
identifications. 

We eventually ignore the efficiency loss due to photo- 
metric selection, and integrate the effects of possible con- 
taminations into a drift of absolute magnitude (see Section 
EIT.3.41 and Eq. [5]). Note that if spectroscopic redshifts of 
host galaxies are acquired and the wrong host is assigned 
to a supernova, this event will likely fail the photometric 
typing cuts and hence will not pollute the Hubble diagram. 

3.6. Rest frame I- band Hubble diagram 

The wide survey provides the possibility of constructing a 
rest frame I-band Hubble diagram out to z ~ 0.9, with 
a statistics of about 7000 events. Encouraging cosmologi- 
cal results were recently obtaine d from rest frame I-b and 
measurements from the ground ([Freedman et al.l I2009T ) . I- 
band distances have a smaller contribution of colour to the 
distance estimate than B-band distances, and are hence 
significantly more robust to systematics related to colour 
modelling and calibration. They also exhibit a smaller scat- 
ter. Measuring distances to the same events using indepen- 
dent measurements constitutes a very appealing test of the 
methods and possible biases. It seems unlikely that collect- 
ing full rest frame I-band light curves from the ground for a 
0(1000) events sample reaching z — 0.9 becomes feasible in 
the next decade. In what follows, we conservatively did not 
include the cosmological constraints that could be obtained 
from this second built-in supernova survey. 

4. Methodology 

4.1. Point source photometry 

In deriving the photometric accuracy of a measurement, we 
assume that the photometry of supernovae is carried-out 
using PSF photometrjj£|. The pixel sampling is accounted 
for, and we checked that the position of the source centre 
within a pixel does not change significantly the signal to 
noise ratio. For most of our supernovae measurements, the 
contribution of the shot noise from the source itself is not 
negligible. 

4.2. Light curve fitter training 

We use the SALT2 model (|Guv et al.l l2007t) to generate 
light curves. In this framework, light curves are described 
by 4 parameters: a date of maximum light (in B-band), an 
overall amplitude of the light curve X$, an X\ parameter 
closely related to stretch factor or decline rate, and a colour 
at B maximum. Since light curve shapes are not derived 
(yet) from explosion models, the light curve models rely 
on a training sample, and then reflect the quality of this 
training sample. Since we should aim at obtaining a super- 
nova sample significantly larger and of better quality than 
before, the sample itself should be used to train the light 
curve fitter. This self-training procedure is only possible for 
light curve fitters which do not provide distance e stimates, 
such as SALT and SALT2 (iGuv et al.ll2005l l2007lV SIFTO 

We assume that the object position is perfectly known: first 
it does not improve the flux variance, and second, we are study- 
ing here a supernova survey where the object is measured at the 
same position in an image series. 
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(|Conlev et alJl2008h . Those do no assume any relation be- 
tween redshift and brightness and only predict light curve 
shapes in observer bands and flux ratios between bands 
or epochs. They leave to a later stage the event param- 
eter combination that will be used as a distance estima- 
tor. Supernova model training causes specific uncertainties: 
model noise due to the finite training sample and the cal- 
ibration uncertainties of the latter, and a specific distance 
uncertainty pattern due to self-training. 

SNe la exhibit a variety of rest frame colours, but dif- 
ferent colours of the same event seem closely related. To 
account for correlated colour variability, one can take guid- 
ance from light extinction (by e.g. dust) and model the 
observed magnitude (for example at time of maximum B 
light) as: 

m bs(X) = mo + mi(X) + c Q(X) (1) 

mo and c vary from event to event, whilst the "naked" flux 
mi and the colour law Ci constitute the model (suppos- 
edly for all events). We can chose the scaling of C; so that 
c = (B — V) max + constant. C'i (called the "colour law" 
in the SALT2 parlance) determines the colour relations, 
e.g. the slopes measured in rest frame colour-colour planes. 
When a similar model is developed for extinction by dust, 
the equivalen t of the Ci function is determined from data. 
For example ICardelli et al.l (1989) propose a one-parameter 
family of extinction laws commonly called "Cardelli laws" . 
For a supernova model, C; can also be extracted from data, 
however up to an arbitr ary constan t which does not affect 
colour-colour position^. iGuv et aT] (|2005l 120071 ) derive C/ 
for SNe la and find it significantly different from Cardelli 
laws. Similarly, the SNe rest fram e colou r relations (e.g. U- 
B vs B-V) found in lConlev et all (p008h are different from 
the ones expected from Cardelli laws. 

The Cardelli laws not only predict colour-colour slopes, 
but they also relate colour variations to total extinction. 
Extinction causes a brighter-bluer relation, characterised by 
i?B=4.1 on average for the Milky Way dust. The brighter- 
bluer relation can again be measured for supernovae (by 
minimising the Hubb le diagram res iduals), and one finds 
a value around 3.2 (|Guv et al.ll2010h . significantly smaller 
than 4.1. We will come back to this point in i J4.3.4l 

Since both SNe colour relations and total to selective 
extinction point away from extinction laws for Milky Way 
dust, we should make provision in the survey design for 
studying both aspects using supernovae, and in particu- 
lar be in a position to measure both colour relations and 
the brighter-bluer correlation without any assumption nor 
prior. It is essential not to restrict the possibilities of SN 
colour variations to those described by the Cardelli laws. 
Even if dust extinction parametrisations adequately de- 
scribed supernova colour relations, the known variations of 
dust properties would still require precise supernova mea- 
surements to define which dust is in action: the fact that 
"regular" dust is unlikely to be the first cause of supernova 

6 The Ci(X) family reported in ICardelli et all (jl989f ) for ex- 
tinction by dust are also obtained from colour-colour slopes. 
The arbitrary constant can however be determined by requiring 
that C;(A) is for A going to infinity. Observations extending to 
~ 5 /im allow the authors to carry out the extrapolation. This is 
not possible in the supernova framework and C;(A) is extracted 
from colour-colour relations up to an additive constant. SALT2 
chooses Ci(\b) = 0. 



colour variations does not really impact on the requirements 
of a supernova survey. The causes of SNe la colour variabil- 
ity are unclear, and collecting as many colours as possible 
for each SN event will likely improve our understanding of 
these matters. 

As noted above, the data scatters around these empir- 
ical modelling of S Ne la colour rela t ions b eyond measure- 
ment uncertainties (|Guv et al.ll2007ll2010t ). We referred to 
this extra noise as colour smearing and used a value of 
(T c = 0.025 mag to set the depth of the observations. This 
finite colour smearing might indicate a fundamental limita- 
tion of extinction-inspired models to describe SNe la colour 
variability. 

4.3. Fisher Matrix 

We want to forecast uncertainties of cosmological param- 
eters expected from various survey scenarios, and various 
hypotheses regarding sources and size of systematic uncer- 
tainties. The Fisher matrix framework is sufficient for our 
purpose. 

Sources of uncertainties can be regarded either as extra 
noise or as extra parameters, and the choice is only a mat- 
ter of convenience, as discussed on a practical example in 
Appendix [X] We settle for adding parameters correspond- 
ing to uncertainty sources, and consider 6 parameter sets: 

1. the cosmological parameters. 

2. the parameters of the SN events themselves. 

3. parameters describing the (rest frame) colours of an av- 
erage supernova. 

4. parameters describing the colour law, i.e. how these 
colours change from event to event. 

5. parameters describing the brighter-slower and brighter- 
bluer relations, and the intrinsic brightness of a super- 
nova. 

6. the photometric zero-points of the light curve measure- 
ments, i.e. how instrumental fluxes are converted to 
physical fluxes. 

We use least squares estimators. Usually, supernova cos- 
mology fits are carried out in two successive steps: the mea- 
sured light curves are first fit to yield the light curve pa- 
rameters, and the cosmology is then fitted using these light 
curve parameters (and their uncertainties). We use here a 
mathematically equivalent approach where a single simulta- 
neous fit considers all parameters at once. The Appendix [Al 
discusses why these two approaches are strictly equivalent. 
We settle here for the simultaneous fit because the same 
events are assumed to be used both to train the light curve 
fitter and to measure cosmology. In particular, photomet- 
ric calibration uncertainties affect cosmology both directly 
and through the light curve fitter training. Within a simul- 
taneous fit, propagating the uncertainties (with correlations 
between events) does not require any particular care. One 
obvious drawback of the simultaneous fit is that considering 
the supernova event parameters on the same footing as all 
the other ones considerably increases the size of the least- 
squares problem and uncomfortably lengthens the require d 
matrix inversion (as experienced in iKim fc Miguel l2006h . 
In Appendix [B] we describe how we take advantage of the 
specific structure of our least squares problem to rephrase 
the linear algebra using only small-sized matrices. 

We now enter into the details of the chosen parametri- 
sations. 
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4.3.1. Dark energy parameters 

We follow the commonly use d equation of state (Eo S ) effec - 
tive parametrisation due to IChevallier fc Polarskil (| 2 lh : 

w(z) = Wq + W a z/(l + z). 

4.3.2. Event parameters 

We use the SALT2 parameters: an overall amplitude, a pa- 
rameter very close to stretch-factor, a date of maximum 
(B-band) light, and a rest frame colour at maximum. 

4.3.3. Average magnitude of a supernova and colour law 

Following the SALT2 model and Eq. [TJ we model the peak 
magnitudes of supernovae as: 

m(A) = m* B + mi (A) + cQ (A) (2) 

where A refers to a rest frame wavelength. The two func- 
tions mi and C; describe respectively the magnitude of an 
average supernova as a function of rest frame wavelength 
and the colour law. We chose, mi (As) = 0, Ci(Xb) — and 
C-i(Xv) — — 1, so that m* B represents the actual magnitude 
in rest frame B band, and c = B — V + constant is the 
colour of the event. Determining mi and C; from the data 
is part of the light curve fitter training. Rather than model- 
ing the functions themselves, we model offsets to the same 
functions extracted from the SALT2 model, by defining: 

mi (A) = m ltS ALT2(X) +£>i(A) 

C l (\) = C l ,SALT2(\)+P2(\) (3) 

where p\ and pi are polynomials with N± and N2 pa- 
rameters to be fitted respectively. The role of these two 
parametrised functions is to correct broadband inaccura- 
cies of an otherwise correct "narrow band" modelling of 
supernovae and their colour variations. These inaccuracies 
are typically expected from photometric calibration errors 
of the training sample, and do not require sharp corrections. 
We settle for 10 parameters for each polynomial, which cov- 
ers spectral resolutions coarser than ~ 20. Our results do 
not degrade significantly with larger values, up to ~ 20 pa- 
rameters, where the Fisher matrix becomes numerically sin- 
gular, because the broadband data does not contain enough 
information at these higher spectral resolutions. 

These two polynomials can be defined either as affecting 
the amplitude of a light curve (at the central wavelength of 
the considered filter) or as multiplicative factors affecting 
the supernova flux before integration in the passband. We 
eventually settle for the second approach because it is more 
realistic. 

4.3.4. Distance estimator parameters 

The simplest way to model the brighter-slower and 
brighter-bluer correlations is via linear coefficients in th e 
distance modulus (see e.g. lTripp|[T99a lAstier et al J 1 20061) : 

H = m* B + aX x - (3c- M (4) 

where X\ may either be the actual X\ parameter from 
SALT2 or any other empirical parameter describing light 
curve width. Ai is the intrinsic B-band magnitude of a su- 
pernova of null X\ and c. Note that a rest frame band dif- 
ferent from B can be chosen, as well as a different colour, as 



in iFreedman et al.l (|2009f ) . All supernova cosmology works 
regard the brighter-slower relation (a) as empirical and 
explicitly or implicitly fit for it. Putting the /3 parame- 
ter on an equal footing is much less common because the 
assumption that colour variations among supernovae are 
caused by dust extinction only readily provides a value 
(or a range of values) for j3 together with a prediction 
for the Ci function above. Many supernova works assume 
that the brighter-bluer relation of SNe, parametrised by 
our j3 parameter, is due to redenning by dust similar to 
Milky Way dust. Although, one should indeed expect that 
some amount of dust extinction is at play, other sources 
might dominate colour variations among events. Fitting for 
/3 constitutes a more general approach than assuming its 
value, especially since th e found val ues depart from Milky 
Way dust red ennin g fe.g.lTripdll998l and ref erences therein 
l Astier et all 120061: IFreedman et all l2009t iKowalski etafl 
120081 lAmanullah et al.ll2010l) . Assuming that colour varia- 
tions of supernovae are due to Milky Way like dust not only 
increases the scatter of distances but also m ay cause arte- 
facts suc h as evidence for a H ubble bubble (|Conlev et al.l 
I2007D . In lAlbrecht et all (|2009h . the difference between the 
expectation that the brighter-bluer relation of SNe follows 
redenning by Milky Way dust and evidence from measure- 
ments that it is not the case is described as a systematic 
uncertainty. We will not follow this route: since distances 
depend on this measurable j3 parameter (or its equivalent in 
a more complex distant estimator), it has to be measured, 
whatever the source of colour variation is. 

In the distance modulus sketched above, a and /3 may 
be either global coefficients or depend on redshift, or de- 
pend on the environment of the supernova (characterised 
for example by host galaxy colours) . We also emulate a pos- 
sible unnoticed evolution of supernovae intrinsic brightness 
(or a smooth redshift dependent distance bias) by allowing 
a redshift dependent intrinsic brightness: 

M(z) = M + e M z (5) 

where A4q and eu are parameters and eu is constrained 
by a Gaussian prior (around a null fiducial value). We will 
consider later several setups for a, j3 and M.. We will use a 
fiducial value of c(ejvr) = 0.01, and we sketch a scheme 
to constrain it from the survey in A5.2I This modelling 
of the "uncertai nty floor" is significantly different from the 
assumptions in iLinder fc Huteren (|2003f ). where distance 
shifts in redshift bins are assumed independent, and the 
uncertainty in a given redshift bin depends on the highest 
redshift reached by the survey. 

4.3.5. Photometric calibration 

Photometric calibration accuracy is naturally a major is- 
sue for supernova surveys. Supernova cosmology "only" re- 
quires relative calibration, in the sense that the heart of 
the method consists in comparing the flux of events across 
redshifts: cosmological results are insensitive to the overall 
flux scale. Supernova fluxes are however to be measured in 
different bands: distant supernova should be measured in 
redder observer bands than nearby supernovae. Since su- 
pernovae fluxes are calibrated against standard stars, the 
calibration uncertainty arises in a first place from our lim- 
ited knowledge of the fluxes of these standards, more pre- 
cisely of the ratio of their fluxes in different bands. A sec- 
ond contribution to flux calibration uncertainty arises from 
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the measurement process itself, i.e. the systematic accu- 
racy of supernova to standards ratio measurement. We will 
examine later the influence of the uncertainty of photomet- 
ric calibration zero-points, where these uncertainties should 
account for both sources. We typically assume that pho- 
tometric zero points are known to 1% including the con- 
version to fluxes. This is conservativ e in view of currently 
obt ained precisions i n the visible (see lRegnault et aTl feOOQ. 
and lGuv et al.ll2010l for the total uncertainty). 

4.4. Cosmologies I priors 

To forecast the measurement precision of supernova sur- 
veys one has to complement the distance measurements 
by other constraints, because distances alone cannot effi- 
ciently separate the various universe densities and the equa- 
tion of state parameters. However, the parameter combina- 
tion probed by distances to moderate redshifts make them 
unique. We will complement the proposed measurements 
with "Planck priors" . For a geometrical probe like distances 
to SNe, these Planck priors essentially consist in a single 
constraint on the geometrical parameters (Om i ^de, w , 
w a ). Supernovae distances yield a second strong constraint 
in the same parameter space, which is not enough to re- 
duce the allowed region in the wo,w a plane. Planck priors 
are more efficient to complement BAO surveys for dark en- 
ergy because CMB constrains a combination of geometrical 
parameters and sets the size of the BAO "standard ruler" . 
As an extra constraint for our supernova survey, we settle 
for the simplest one: flatness. In using CMB priors for dark 
energy, one should pay attention to not making use of dark 
energy information available through the ISW effect, be- 
cause this information available on large scales might not 
be reliably extracted. This is achieved by exactly enforc- 
ing the "geometrical degen eracy" (see e.g. lAlbrecht et al.l 
2009; lMukheriee et al.ll200ll . Within this framework, using 
the full Fisher matrix or a single one-dimensional geomet- 
rical constraint makes little difference. Our "Planck priors" 
hence reduce to a measurement of the shift parameter R 
(R = \/£LMHnr(z C M B)) to 0.32% relative accuracy (see 
iMukheriee et al.l l2008. Table 1), together with flatness. 

5. Results 

We simulate a fiducial flat ACDM universe with Om = 0.27. 
The equation of state is parametrised as w(z) = w + 
w a z/(l + z). We assume that the scatter of the Hubble 
diagram for perfect measurements exactly matching the av- 
erage supernova model (i.e without colour smearing) is 0.12 
mag nitudes. Current estimates of this quantity are below 
0.10 (|Guv et al.ll2010D . but since it can only be obtained by 
subtraction of identified uncertainties that might have been 
inadvertently inflated, we choose to stand on the safe side. 
We assume that the colour smearing (see i j3 . 1 [) causes peak 
SN magnitudes in a single ban d to scatter aroun d the model 
by a quantity a c = 0.025 (see lGuv et al.|[2010l) unless oth- 
erwise specified. a c — 0.01 was assumed in iKim fc Miguel 
(1200611 . 

We restrict the rest frame central wavelength of the 
bands entering the fit to [3800-7000] A, which leaves 3 to 4 
bands per event. Enlarging this range formally improves the 
performance but breaks the requirement that similar rest 
frame ranges are used to derive distances at all redshifts. 



In a real survey, the whole information would of course be 
used, in particular to study the supernovae colours. 

We will now study the survey performance regard- 
ing constraints of the equation of state. Following 
lAlbrecht et al.l (|2006l ). we define the pivot redshift z p as the 
one where the EoS uncertainty is minimal, and w p = w(z p ). 
As performance indicators, we report a(w p ), the uncer- 
tainty on of the EoS evolution a(w a ), and the reciprocal 
of their producfl often used as a figure of merit (FoM). 
The quantity o~(w p ) can be regarded as the ability of the 
project to challenge the cosmological constant paradigm. 

In Table [5J we turn various uncertainty sources on and 
off labelled from A to Z. The setup A only considers pho- 
tometry Poisson noise and Hubble diagram scatter, and the 
baseline scenario Z adds all other uncertainty sources dis- 
cussed above. Some combinations of uncertainty sources are 
provided in between, in order to identify major uncertainty 
drivers. Lines B to F display the effect of one source at a 
time : the single source that mostly degrades the perfor- 
mance is the intrinsic brightness drift with redshift cr(eM), 
followed by the photometric calibration (z p ). In lines G to 
M, both intrinsic drift with redshift (a(eM) = 0.01) and 
colour smearing (er c = 0.025) are allowed, and we add 
other sources. We note that adding the supernova model 
fit alone (lines H,I,J compared to G) does not significantly 
degrade the performance and that adding in the calibration 
uncertainty (K compared to G) has a sizable effect. But the 
most dramatic effect is the combination of supernova model 
training and calibration uncertainty (lines H,K and L), be- 
cause the Ni parameters of pi (Eq. [3]) and the zero-points 
play the same role, the first set in the supernova frame and 
the latter in the observer frame. This enlights the role of cal- 
ibration uncertainty through the light curve fitter training, 
and is s upported by the findings of current ground-based 
surveys (|Guv et alj|20Tol) . Turning on or off the colour law 
fitting (N 2 , comparing lines G and I, or L and Z) has es- 
sentially no impact on the performance : there is indeed no 
benefit to rely on an assumption (such as a Cardelli law) 
for this part of the model. 

One might be surprised that calibration uncertainties do 
not ruin the proposed survey since current ground-based su- 
pernova surveys face calibration-induced uncertainties com- 
parable to statistics with "on ly" a few hundred events 
(jGuv et all [20101 IConlev et al.l(2010D . There are two key 
differences: first, these survey have to face cross-calibration 
issues with the nearby sample; second, the observed rest 
frame region varies with redshift, and the comparison of 
events across redshifts then heavily relies on the supernova 
model and inherits its uncertainties. This illustrates our 
requirement that all events be observed in the same rest 
frame range. 

Comparing colour distributions along redshift consti- 
tutes a sensitive handle on astrophysics conditions evolv- 
ing with time. This test is independent of the supernova 
model training, in the sense that the training does not aim 
at matching these distributions. The photometric calibra- 
tion accuracy directly limits the comparison of colours at 
different redshifts: with zero-points defined at 0.01, colours 
can only be compared to ~ 0.015. Given the natural colour 
spread of ~ 0.1, this is a test at a/7 level, which does not 
benefit from more than ~ 50 events in a redshift slice. For 



7 This quantity is equal to Det(C ov(wo,w a )) and scales 
as the figures of merit inspired from I Albrecht et al.l (2006). 
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label 


Nt 


N 2 


cr(zp) 


0" c 


er(eM) 


cr(w a ) 


z p 


a(w p ) 


\a(w a )a(wp)] 1 


A 

















0.23 


0.41 


0.020 


213.9 


B 


10 














0.23 


0.41 


0.020 


213.0 


C 





10 











0.23 


0.41 


0.020 


213.9 


D 








0.01 








0.25 


0.40 


0.021 


193.0 


E 











0.025 





0.25 


0.41 


0.020 


198.1 


F 














0.01 


0.28 


0.35 


0.024 


152.1 


G 











0.025 


0.01 


0.29 


0.36 


0.024 


141.7 


H 


10 








0.025 


0.01 


0.29 


0.36 


0.024 


140.5 


I 





10 





0.025 


0.01 


0.29 


0.36 


0.024 


141.7 


J 


10 


10 





0.025 


0.01 


0.29 


0.36 


0.024 


140.5 


K 








0.01 


0.025 


0.01 


0.32 


0.37 


0.027 


117.0 


L 


10 





0.01 


0.025 


0.01 


0.42 


0.32 


0.031 


76.9 


M 





10 


0.01 


0.025 


0.01 


0.32 


0.37 


0.027 


116.9 


N 


10 


10 


0.01 


0.025 





0.40 


0.33 


0.030 


83.7 


Z 


10 


10 


0.01 


0.025 


0.01 


0.42 


0.32 


0.031 


76.8 



Table 5. Dark energy constraints when considering specific subsets of uncertainties. Ni and N2 refer to the number of 
parameters describing the supernova model (Eq. [3]), u(z p ) is the photometric calibration uncertainty for each band, <j c 
is the scatter of light curve amplitude around the model, and a(eu) is the prior applied to the evolution parameter cm 
defined in Eq. [SJ z p denotes the (pivot) redshift which minimises the equation of state uncertainty, and this minimal 
uncertainty is labelled a(w p )). The last column is the reciprocal of the area of the marginalised uncertainty ellipse 
(A% 2 = 1) in the (wo,w a ) plane (and also in the (w p ,w a ) plane, see lAlbrecht et al1 [2~006. p. 97). 



supernovae, the sensitivity of evolution tests constitutes an- 
other strong incentive to improve the calibration precision 
to a few per mil level. 

5.1. Variations of performance when altering parameters of 
the baseline. 



Alteration 


er(w p ) 


[o-(iy a )o-(iy p )] _1 


none 


0.031 


76.8 


Ni = 15 


0.031 


75.4 


N 2 = 15 


0.031 


76.8 


(Tint = 0.10 


0.029 


86.7 


<r{e M ) = 0.02 


0.033 


64.6 


a{zp) = 0.005 


0.028 


99.3 


cr c = 0.015 


0.029 


86.3 


half statistics 


0.037 


50.6 


double statistics 


0.025 


118.1 


only wide 


0.044 


39.2 


only deep 


0.046 


31.9 


wide x 0.82, deep x 1.25 


0.031 


74.7 


wide x 1.18, deep x 0.75 


0.031 


78.0 



Table 6. Effect of altering survey parameters of the base- 
line (defined as line Z of Table [5]). 



Table [5] illustrates the effect of parameters that deter- 
mine the cosmological performance, and varies them one 
at a time. We can check that increasing the spectral res- 
olution of the supernova model (N\ and N2) has almost 
no effect. If an intrinsic resolution <Ji n t better than 0.10 is 
confirmed, this yields a ~ 10 % improvement of the FoM. 
Doubling the evolution systematics (J {cm) degrades it by 
16 %. Improving the zero-points accuracy by a factor of 2 
improves it by 30 %. The joint effect of calibration uncer- 
tainty and distance biases on a(w p ) is displayed in Fig. [5] 
One might also note that reducing statistics significantly 



alters the performance. Conversely, doubling the statistics 
improves the FoM by 50 %. We should probably stress here 
that the wide survey should not be regarded as doable from 
the ground, because supernova distances at z ~ 1 make use 
of y and J bands. We finally vary the wide and deep survey 
allocations within a constant overall observing time, and 
note that this is not a key parameter. Reducing the deep 
survey marginally improves the FoM but at the expense of 
reducing the high redshift statistics, which is the most pre- 
cious to tackle evolution issues. One might also regard the 
rest frame I-band Hubble diagram as an extension of the 
wide that would improve the constraints. If, following the 
SNAP approach, one includes 1000 nearby supernovae at 
z = 0.05 (measured from the ground), assuming that cross- 
calibration uncertainties are not worse than assumed here, 
and ignoring potential bias issues, the FoM reaches 110. 

One should seriously consider the possibility that there 
are subclass es of SNe la depending on environment as 
sugge sted by iMannucci et all (|2006ft . ISullivan et all (|2006l 
12010ft propose some observational evidences for different 
average properties of SNe occurring in passive and star- 
forming galaxies. It is not yet clear if these different envi- 
ronments produce different supernovae, or if these differ- 
ent environment sample differently the same parent pop- 
ulation. We will consider however here the most dramatic 
case, where two environments produce two different event 
species described by different parameter sets, namely dif- 
ferent a, (3, M (see equation 2]) and different supernova 
models (different p\ and pi functions, see equation [3]) . We 
consider that the host galaxy colours allow one to assign 
each event in one category. In the case of the admixture 
not evolving with redshift and categories having the same 
photometric quality, the variance of the cosmological esti- 
mators is mathematically the same as for a single species 
scenario, as shown in Appendix [Cj For a more realistic sce- 
nario, we varied the admixture with redshift from 30/70% 
at z = to 70/30% at z = 1.5, and variances of the cosmo- 
logical parameters do not increase by more than 1%. 
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0.03 



a(e ) 

Fig. 2. Contour levels of a{w p ) as a function of calibration 
uncertainty ozp (equal for all filters), and the distance evo- 
lution uncertainty cr(eA/) (defined in Eq. [5]). 



So, if it turns out that SNe la consist in an admixture 
of species that can be tagged using host galaxy colours or 
properties of the supernova light curve, the proposed strat- 
egy retains its performance. Note that even if events are 
assigned the wrong category, the cosmological parameters 
remain unbiased as long as the fraction of interlopers does 
not evolve with redshift. We also considered the possible 
evolution with redshift of a and j3 parameters (defined in 
equation [4]) by fitting those in redshift bins of 0.1, and vari- 
ances are changed at the percent level. Even considering 
simultaneously event sub-classes and redshift bins results 
in minute degradation of performances. 

Finally, we replaced the evolution of A4(z) defined in 
Eq.0by: 



M(z) = M + e M 



[dLogd^/ Owq] (z) 
[dLogd L /dw ](z = 1) 



where the denominator has been chosen so that a (cm) rep- 
resents the uncertainty at z — 1 as in Eq. [5] This functional 
dependence makes ejvf degenerate with Wq. The choice of a 
sys tematic mimicking the dep endence on wq was inspired 
by lAmara fc~ R cfrcgier (2008|) where, in the weak shear 
framework, it is regarded as a worst case. However, in our 
case, for the same strength of the prior, we find that the 
FoM are larger (or equal) than with the linear model of Eq. 

El 



5.2. Constraining a redshift dependent bias of distances 

As discussed in the previous paragraph, cosmological con- 
straints from supernovae depend on our ability to constrain 
the systematic drift of distances with redshift. We investi- 
gate in this section a few handles that might be used to 
bound the effect or correct for it. 

The proposed survey provides accurate colours of host 
galaxies extending to the observer NIR. This enables one 



to estimate host galaxy properties by comparing measure- 
ments to synthetic galaxy spectra, and analyse if derived 
distances (or other supernovae properties) depend on the 
host galaxy propert ies. This is the strategy followed in 
ISullivan et al.l (|2010[) . where evidence for a dependence of 
distances with host galaxy stellar mass is presented, and 
found harmless for cosmology. With the supernova statis- 
tics we are considering here, this approach becomes even 
stronger because it can be applied within a modest red- 
shift range where photometric calibration issues affect all 
supernovae and host galaxies in the same way. 

The metallicity of exploding white dwarves is ex- 
pected to increase (at least on average) with cosmic time. 
Metallicity certainly infl uences the amount of 56 Ni syn- 
thetised in the explosion (|Timmes et al.ll2~003t ) , but we have 
no hint yet that standard distance estimators do not cor- 
rect for this evolution. One might even argue that the range 
of environments found at, e.g., small redshifts efficiently 
"trains" distance estimators so that they remain unbiased 
as redshift varies. 

However, in order to bound a possible redshift depen- 
dent distance bias we propose to use the near UV flux vari- 
ations that explosion mode l s correlate to metallicity vari- 
ation s (jHoeflich et al.lll998t ILentz et al.ll2000t ISauer etall 
120081 ). If the models do not necessarily agree on the size 
of the effects, they define the 2500-4000 A spectral region 
as sensitive to admixtures in the progenitor material of 
other elements than Carbon and Oxygen (the ones of the 
baseline scenario). SNe la exhibit a colour diversity, but 
different colours appear to be tight l y connected (see e. g 
lAstier et al1l2006t IConlev et al.ll2005 iFolatelli et al.ll2010D . 
Defining U* as the 2500-3200A region, the rest frame 
combination (U* — B) — 5A(B — V) is, from observations, 
the smallest scat ter combination of the (U*,B,V) triplet 
(|Guv et al.ll2010l) ( up to a mult i plicat ive constant). Using 
synthetic data from iLentz et al.l (|2000t) , we can check that 
this three-band combination is both sensitive to metallic- 
ity, and nicely correlates to distance biases due to evolving 
metallicities. Quantitatively, when varying metallicity, the 
distance shift <5/i varies as 0.1[(U* — B) — 5.4(S — V)], which 
allows one to constrain 5fi at the 0.01 level in the presence 
of calibration uncertainties of 0.01 per band. The U* band 
is observable in the proposed survey beyond z = 0.6, and 
lies below the supernova model spectral region, and is then 
ignored in the training. 

Beyond the decline rate paradigm, the early phases 
of light curve are expected to encode metallicity (see 
iHoefhch et aDll998t ). Given the envisaged statistics, minute 
departures from the average light curve shape can be de- 
tected and correlated with other observables. One key qual- 
ity of light curve shape indicators is that they remain un- 
affected by calibration issues. 



6. Discussion 

6.1. Comparison with weak shear performance 

The ability of space based measurements of the weak 
shear to constrain the cosmological model, and in par- 
ticular its dark en e rgy s ector was studied in detail in 
lAmara &; Refr egicr (2008) (and references therein). This 
work su pports the st rategy developed for the EUCLID 
project ( Laureiisl l2009[h and identifies the measurement of 
the shear (intimately related to the measurement of second 
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moments of galaxies) as the key systematic. The shear mea- 
surement uncertainty translates to a systematic uncertainty 
floor of the shear angular pow er spectrum C s ( vs . Current 
shear measurement techniques (Bridle et al. 2010) achieve 
a systematic u ncertainty that would limi t a(w p ) to about 
0.05 (Fig 11 of lAmara fc Refregierl (|2008T) interpolated for 
C„ vs ~ 2 10~ 6 from the best algorithm in Table 5 of 
I Bridle et al.| [2010') . One should note that this performance 
is obtained on simulations where in particular, the PSF of 
the imaging system used to measure the shear is perfectly 
known. The EUCLID p roject requires that C s e ys < 10" 7 be 
reached (|Laureiisl|2009 ). leading to <j{w p ) ~ 0.02. 

Using the current performance of ground-based dis- 
tance measurements to supernovae, we conservatively de- 
rive an EoS constrain <t(u> p ) ~ 0.03 from a space-based 
survey. This is in a position to really complement EoS con- 
straints from shear correlations. These two approaches are 
not only complementary because their redshift dependent 
biases are unrelated. On the one hand, distance tests are 
independent of growth rate, and have to be complemented 
(by e.g. Planck priors) in order to constrain the EoS, and 
on the other hand, shear correlation tomography can au- 
tonomously constrain the EoS by assuming a given rela- 
tion (from e.g. General Relativity) between distances and 
growth rate. 



6.2. Summary 

We have proposed a two-cone supernova survey conducted 
from space with a modified EUCLID setup : we assumed 
that both the visible and the NIR imagers are equipped 
with swappable filters. We find that it is possible to accu- 
rately measure more than 10 4 supernovae at 0.15 < z < 
1.55 in 18 months of survey. The photometric accuracy 
is tailored to match the measured intrinsic variability of 
colour relations of supernovae at the highest redshift of the 
survey. All events are measured in the 7 instrumental bands, 
and the BVR rest frame bands are covered at all redshifts. 
Our analysis of supernova distances relies on conservative 
assumptions and the current know how. It integrates many 
nuisance effects, such as the light curve fitter noise together 
with the impact of a conservative photometric calibration 
uncertainty both directly on cosmology and through the 
light curve fitter training. Our approach ensures that the 
interplay of identified uncertainties is properly accounted 
for. The proposed observing strategy also collects the data 
to build a rest frame I-band Hubble diagram to z ~ 0.9, 
with ~ 7000 events. 

Our results are encouraging in the sense that including 
these realistic nuisance effects, competitive constraints of 
the dark energy equation of state can be obtained, when 
using a simple geometrical Planck prior: within a two- 
parameter dark energy model, the EoS can be constrained 
to -0.03 at z ~ 0.3. 

Acknowledgements. We are grateful to P eter Nugent f or pro viding us 
with the simulated spectra discussed in [Lcritz et al. (200TJ), and to 
Andrew Howell for fruitful exchanges. 



Appendix A: Simultaneous versus successive least 
squares fits 

We compare here two approaches for the least-squares fits 
in the supernova framework. In the traditional approach 



one first fits separately the event light curves to extract the 
event parameters, and then fits the cosmology to distances 
constructed from these event parameters. In the less con- 
ventional approach we present here, light curves and cos- 
mology are fitted simultaneously, by minimising the sum: 



(A.l) 



with 



E 



(ji(Pc,Zi)-h*6i-M) 3 



XLC,; 



E 

fc=l,JVi 



(fik - <f>(6i,t ik )f 



(A.2) 
(A.3) 



where, i indexes SN events, k = l..N{ indexes the mea- 
surements fik of a an event, <j> is the supernova model, 
6i are the parameters of event i (9 denotes the ensemble 
of event parameters), p c are the cosmological parameters, 
/J,(Pc, z) is the distance modulus at redshift z, h^Oi is the 
measured distance modulus of event i (a linear combina- 
tion of event parameters) , and M. is a combination of the 
intrinsic magnitude of a supernova and Hq. Oi n t is the in- 
trinsic dispersion of a supernova (defined as the "observed" 
Hubble diagram scatter for ideal measurements), and <Jik 
refers to measurement uncertainties. The first \ 2 term fits 
cosmology, the second fits light curves, and both terms are 
related through 9i parameters. 

If one minimises separately the light curve parts with 
respect to 9i , x\c i mav be approximately re- written: 



Xlc,i = Vi - %) T Wl(9i - 9[) 



(A.4) 



where 9[ minimises Xlc i an< ^ W/ ^ s second derivative 
matrix w.r.t 9i. The approximation would be exact if the 
supernova model were linear with respect to its parameters. 
The approximation holds for error propagation at first or- 
der, and will be used only for this purpose. In order to 
compare the simultaneous fit of p c and all 9i with the two- 
step process, we will compute x 2 (Pc, 9) where 9 denotes the 
9 value that minimises the global \ 2 (|A. 1|) for a given value 
of p c . A tedious calculation yields: 



X 2 (Pc 



N srI . 

E 

i=l 



(fx(p c ,Zi) - h 



To'i 



M) 2 



h T W'h 



(A.5) 



which corresponds to the standard cosmology fit that one 
would perform after a separate fit of all events to yield 
9\ and W[. So, fitting all events and the cosmological pa- 
rameters simultaneously (using expressions IA.2I and IA.3I) 
yields exactly the same cosmological estimators as the "tra- 
ditional" two-step method. Note that the expressions IA.2I 
and IA.5I have in particular different denominators, indicat- 
ing that the simultaneous fit is a convenient way of prop- 
agating uncertainties of other analysis steps. This result is 
very general, and should be no surprise to Kalman filter 
practitioners. 
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Appendix B: The local-global technique 

For a supernova sample of several thousand events, the 
number of parameters in expression IA.1I becomes uncom- 
fortably large (even for simulations), at least when it comes 
to inverting the \ 2 weight matrix. However, since no y 2 
term involves the parameters of two different events, this 
weight matrix is sparse. We write it block-wise: 



W = 



A] 



A x 
W x 


••• Ai 



••• A N \ 






'■. 












'-. 

W N J 



\A T N 

The first row and column refer to p c and the following ones 
to the successive 9i . It is easy to check that the correspond- 
ing blocks of the inverse matrix read: 



N 



c p = (w p -Y,A,w^Ajy l 

i=l 

C Pi i = —CpA^W^ 

Cij = S^W' 1 + W^AjCpA.W- 1 
The solution of the linear system: 



W 



(p\ 

01 


-R 


\0n) 


\B N J 



can be written in the efficient form: 

N 

p = Cp{B p -Y J A l Wr 1 B i ) 

i=l 

e i = W i r 1 (B i -Ajp) 

which only requires the C p block from the inverse. Note 
that the computing time scales as the number of events N, 
rather than N 3 if some standard factorisation of W were 
used. In our case, the C p matrix is obtained in matter of 
seconds with 0(10000) supernova events. 

In order to fit in this scheme, least squares problems 
should exhibit "local" (8i) and "global" (p) parameters. 
By definition, the local parameter sets 0i are not con- 
nected to each other through any y 2 term, and the single 
global parameter set p gathers all the remaining ones. In 
the cases discussed in this work, the split is fairly obvious: 
local parameters refer to supernova event parameters and 
the global set p contains everything else (cosmology, zero- 
points, ...). This block solution scheme of our least squares 
mini misation was propo sed for a totally different purpose 



minimisati on was propos< 
in lReenault et all ([2009). 



Appendix C: Several event species 

We prove here why splitting the supernova sample into sev- 
eral event species does not degrade the cosmological pa- 
rameters, if these species have the same redshift distribu- 
tion and photometric accuracy. We formally split the global 



parameters (all parameters but the supernova events ones), 
into the common ones (labelled by c) which typically gather 
cosmology and zero-points, and the parameters specific to 
each species (labelled by 1 and 2). The later gather the su- 
pernova model and the distance estimator parameters. The 
global parameter space addresses {g c , 31,52), and we wish 
to compute the covariance matrix of g c (which contains the 
cosmology) when fitting both samples simultaneously. We 
add the inverses of both samples covariance matrices, ex- 
pressed in the proper subspaces: 



W 1+2 = W 1 + W 2 = 



x 1 




(w c2 







w 1 


9* 





















w 2 J 



XT 




'W cl +W c2 Xi AY 
ATf Wi 

xl w 2/ 



Using Eq. IB. II the covariance matrix of the common pa- 
rameters reads: 



(B.i) Cov(g c ) = [(W cl + W c2 ) - X^Wr 1 X 1 - XfW^Xi] 1 

If both species have the same properties (redshift distri- 
bution, photometric quality), we have W\ = kW 2 and 
X\ = kX 2 , where k stands for the population ratio. So, 

Cov(g c ) = [(W cl +W c2 )- 

(X l + X 2 ) T {W l + w 2 )-\x l + X 2 )] ~* 

This expression describes as well all events belonging to a 
single species, which completes the proof. As expected, the 
argument breaks down if there are not enough events to 
constrain g\ or g 2 , i.e. if W\ or W 2 is singular. One may 
also note that the argument applies as well to any number 
of event species. 
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